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Abstract. 

y-ray astronomy has produced for several years now sky maps for low photon statistics, non-negligible background and 
comparatively poor angular resolution. Quantifying the significance of spatial features remains difficult. Besides, spectrum 
extraction requires regions with large statistics while maps in energy bands allow only qualitative interpretation. 
The two main competing mechanisms in the VHE domain are the Inverse-Compton emission from accelerated electrons 
radiating through synchrotron in the X-ray domain and the interactions between accelerated hadrons and the surrounding 
medium, leading to the production and subsequent decay of ;r" mesons. The spectrum of the VHE emission from leptons is 
predicted to steepen with increasing distance from the acceleration zone, owing to synchrotron losses (i.e. cooled population). 
It would remain approximately constant for hadrons. 

Ideally, spectro-imaging analysis would have the same spatial scale in the TeV and X-ray domains, to distinguish the local 
emission mechanisms. More realistically, we investigate here the possibility of improving upon the currently published HESS 
results by using more sophisticated tools. 
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1. DESCRIPTION OF METHODS 

The most simple procedure to produce hardness maps 
rehes on bins of fixed size, defined independently of 
the source morphology. The bin size is a compromise 
between increasing the number of excess events per bin 
and keeping small-scale structures. This guarantees the 
absence of any internal artefact but is not adapted in case 
of large flux variations. 
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For 7-like (on) and background (off) events in two sepa- 
rate high and low energy bands, these formula define the 
excess E iQl], the contrast C (used here instead of the im- 
mediate ratio of excesses) between the bands and its sig- 
nificance Sc, combining the significance in each band, 
given by the Li&Ma formula [|2j]. The selection of the 
energy bands, starting at the energy threshold for the ob- 
servation, aims at balanced counts of excess events. 
Adaptive binning 

In the Weighted Voronoi Tessellation (WVT) adaptive 
binning algorithm jsl], the core procedure assigns pix- 
els to non-overlapping bins, whose shape is maintained 



close to circular From a given position Pstart, neighbour- 
ing pixels are accreted until the requested significance 
or maximum size is reached. The assignment of pixels 
is repeated to reduce the spread of bin significance and 
shapes. Identified point sources are treated here as fixed 
bins, to avoid contamination by the background as well 
as smearing out of the source. To make the whole al- 
gorithm independent of the start position, we loop over 
random Pstart, keeping the average of the re-binned maps. 
This procedure is illustrated in Fig.|2] The target signifi- 
cance is 7 (7, with an upper limit on bin radius of 0.12°. 
Adaptive smoothing 

Adaptive smoothing methods, like ASMOOTH 101 algo- 
rithm (CSMOOTH in the CIAO package iH), adjust the 
scale of a smoothing kernel to maintain a constant sig- 
nificance. For a given scale, the image is smoothed and 
pixels reaching the desired significance are added to the 
output image and removed from the input one. Remain- 
ing pixels are treated with progressively larger scales. To 
build the contrast map, the two energy bands are treated 
separately, see 01 for details, with a 9 (7 target here. 
Photon index 

From the contrast, we estimate the photon index and flux 
normalization, assuming a power-law spectrum. A hard- 
ness ratio derived from excess counts instead of fluxes is 
affected by the instrument response. To correct for this, 
the effective area of the instrument is integrated in the 
energy bands, accounting for zenith angles and offsets 
from the observations. 



FIGURE 1. Applying different methods to tlie same count map. From left to right, fixed-scale smoothing, re-binning, AS- 
MOOTH. The contours are for significance, from 3 to 7 <7 or 5 to 9(7 for ASMOOTH). 




FIGURE 2. The WVT algorithm. From left to right, the original WVT, the adapted WVT and the associated average bin sizes. 
Significance contours from 3 to 7 a. 
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FIGURE 3. Slices for model I, with the geometry of Figure \Qeft. In black is shown the model. Green, blue and pink are 
respectively the fixed binning (0.04° integration radius), WVT and ASMOOTH algorithms, with 1 a error bands. 



The acceptance profiles are built for the same energy 
bands from regions in the sky devoid of 7-ray source. For 
each position on the map, the spectral power-law index is 
varied until the hardness ratio agrees with the data. The 
flux is then scaled to the excess. 



Disclaimer 

Several important issues are not investigated here, no- 
tably the variation of the PSF with energy (smaller at 
higher energy) and its effect on the actual spatial resolu- 
tion, nor the inclusion of trial factors in the significance. 
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FIGURE 4. Slices for model II, following the conventions of Fig[3] 



2, APPLICATIONS 

Simulating a shell-type supernova remnant 

RX J1713-3947 is a large shell-type SNR with bright 
non-thermal X-ray emission |0, [31 ■ Recent XMM lit] 
and Nanten ^ observations suggest interactions with 
molecular clouds, at ~ 1.3 kpc. TeV spectro-imaging 
was obtained by HESS ifiol [Till , showing good correla- 
tion with the X-ray morphology. The TeV spectral in- 
dex seems constant across the remnant, at odds with 
X-rays The nature of the TeV emission is still de- 
bated IHtil [US- 

We investigate the visibility by HESS of two qualita- 
tive models: one (model I) with a spectral photon index 
decreasing from the shock front (-2) to the center (-2.4), 
the other (model II) with a constant index (-2.1). The flux 
distribution follows a shell-type profile (spherical geom- 
etry) with a gradient from South-East to North- West. 50 
simulations, mimicking the current HESS statistics in the 
energy bands 0.45 - 1 TeV and 1.2-30 TeV , were pro- 
duced to estimate the fluctuations. 
The contrast profiles, Fig.|3]and|4] are properly recovered 
by the three methods. The significance values are also as 
expected: varying with the flux level for the fixed binning 
and stable near the requested level everywhere within the 
emission region for the other methods. For the same rea- 
son, small scale details are lost, as shown for model I by 
the biased contrast given by ASMOOTH and WVT at the 
center of the shell. This suggests that fixed binning may 
have more potential for the discovery of sources or emis- 
sion features, to be confirmed by additional observation. 
Overall, from the errors on the contrast, the variation for 
model I is only a 2-3 (7 effect on average. Discriminating 
between the two models included here would probably 
be difficult on a real source. 



Pulsar wind nebula HESS J1825-137 

HESS J 1825- 137, detected in the HESS Galactic Sur- 
vey iH], is likely associated with the X-ray PWN G18.0- 
0.7 and the very energetic (3x10^^ erg/s) pulsar PSR 
J1826-1334 inil. The VHE emission is offset to the 
South / South- West of the pulsar, y-ray spectral indexes 
m the PWN iH soften with increasing distance to the 
pulsar. 

From the energy bands 0.27 - 0.9 TeV and 1.1 - 30 TeV, 
we produced maps of the photon index. Fig. |5] and the 
flux density, Fig.|6l The estimation of errors are prelimi- 
nary. 

Our results on the PWN are compatible with the HESS 
publications: the evolution of the photon index is clearly 
seen, with a similar range of values, and so is the flux 
normalization profile. The superposed significance con- 
tours indicate that the external part of the profile is not 
reliable. 

The neighbouring binary LS 5039, a point source for 
HESS, was treated as a fixed bin of radius 0.1 °. The in- 
tegrated flux and photon index are here again very close 
to the published |[l9ll phase-averaged numbers. 

3. CONCLUSION 

The methods discussed here give results compatible with 
HESS publication. They may help exploring the mor- 
phology of extended VHE sources, although as illus- 
trated in a simulation, the uncertainties might prove too 
large for detailed studies. 
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FIGURE 5. Photon index from the WVT maps. The significance contours (black) in the map are for 3 to 7 cr). The values in the 
histogram correspond to the concentric slices (white) in the map, centered on the pulsar. 
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FIGURE 6 

defined similarly to Fig.[5l 



Photon flux density from the WVT maps. The map is scaled in ph. cm ~s 'deg ~. The contours and histogram are 
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